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Abstract. Using 2D nonlinear simulations, we study the 
generation and nonlinear evolution of spiral structure in 
a star-forming multi-component gravitating disk. We con- 
firm in agreement with previous studies the destabilizing 
role of a cold gaseous component and extend this con- 
clusion for multi-component star-forming disks exchang- 
ing mass and momentum between its components. We 
show that the spiral structure growing on a non-stationary 
multi-phase background reaches its saturation in a similar 
manner like the one-component disks. The spiral structure 
survives even if most of the gas is transformed into stellar 
remnants of larger velocity dispersion. 
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1. Introduction 

The regular spiral arms observed in nearby disk galaxies 
are a manifestation of global spiral modes developed in 
galactic disks. Numerous linear and nonlinear studies have 
demonstrated that self-gravity plays a major role in the 
generation of global spiral modes (see, e.g. Binney and 
Tremaine [T987| ). 

Growth of spirals on early stages of galactic evolution, 
however, occurs on a rapidly changing background. The 
precursors of present systems look bluer than old galaxies 
giving evidence for an enhanced star formation rate in 
early epochs. The colors and irregular structure seen in 
intermediate red-shifted galaxies suggest that they are a 
population of galaxies with massive star formation rapidly 
consuming gas in the galactic disks. These galaxies are the 
systems where star-formation might considerably change 
the physics of spiral dynamics. 

In this paper we address the question how the evolu- 
tion of the large-scale spiral perturbations is affected by 
the processes of mass and momentum exchange in a star 
forming disk. Kato ( 1972 1974 ) discussed the effects of 
star-gas mass exchange in a linearly growing spiral mode. 
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Using the WKB-approximation he found, that under some 
conditions the mass exchange between stars and gas can 
excite density waves. We extend this study by following 
the dynamics of the nonlinear hydrodynamical equations 
for three interacting components: the massive stars, the 
low mass stars (or stellar remnants) and the gas. We per- 
form a series of 2D multi-phase hydrodynamical simula- 
tions of the dynamics of unstable multi-component disks, 
and compare our results with the morphological proper- 
ties of spirals, growing in one-component self-gravitating 
disks. As a particular realization of the mass transforma- 
tions in star-forming systems we choose the description of 
Koppen et al. ( |1995| ) which is a subset of the more de- 
tailed " chemo-dynamical" models developed by Theis et 
al. ( |1992D , Hensle r et al. ( |1993D , Samland ( |1994| ) and Sam- 
land et al. ( |1997| ). Apart from previous papers investigat- 
ing the interaction network of a multi-component system 
for spherical (Theis et al. 1992) or axisymmetric galax- 



ies (Samland et al. 1997) we deal here with the evolution 
of thin, but non-axisymmetric systems. This allows us to 
study the influence of the interactions between the com- 
ponents on the evolution of spiral structure. 

2. Basic Equations for a Three-Component Star 
Forming Disk 

The physical processes in star forming galactic disks are 
too complicated for a detailed analysis from first princi- 
ples. Among numerous dynamical and interchange pro- 
cesses one has to choose those which have a major in- 
fluence on the evolution of the system. We perform our 
analysis using a simplified model which splits the disk 
into three components regulated by time-dependent mass 
transformations. We take into account a gas component 
containing all phases of the interstellar medium, and two 
stellar components, the massive and the low mass stars. 
The less massive stars are assumed to have no influence 
on the interstellar medium, and are simply accumulated 
as remnants, whereas massive stars are assumed to be re- 
sponsible for gas heating. 

The chain of mass transformations within the model 
includes spontaneous star formation and ejection of the 
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stellar mass back into the gas phase. For the description 
of the star formation rate we use the approach developed 
by Koppen et al. (1995) which is a basic skeleton of the 
interaction scheme used in chemo-dynamical models (e.g. 
Theis et al. |1992| , Samland et al. |1997D . It describes the 
spontaneous stellar birth as a power law function of the 
gas density with an efficiency depending on the tempera- 
ture of the gas. The latter just has to guarantee that the 
efficiency drops down with increasing temperature. In this 
scenario, the chain of mass transformation processes can 
be described by the following set of equations: 



dt 



-'^ig,T) + r,s/T 



J = C*(,,T)-./r 
-^il-C)^ig,T) + {l-rj)s/T 



(1) 



(2) 



(3) 



where the star formation rate ^'(5, T) is given by a power 
law dependence on the gas density g and a function of the 
gas temperature T 



^{g,T)^Cnf{T)g\ 



(4) 



s denotes the density of massive stars, whereas r corre- 
sponds to the low-mass stars and the stellar remnants. The 
parameters ( and rj determine the fraction of newly born 
massive stars, and the fraction of gas ejected by massive 
stars into the ISM, respectively. The parameter r corre- 
sponds to the mean stellar lifetime of a massive star. 

In general, for the determination of the "efficiency" 
factor /(T) which depends on the gas temperature, one 
has to in tegrate the energy balance equation. Koppen et 
al. (1995) noticed, however, that the model can be simpli- 
fied, if the time-scales for heating and cooling of the gas 
are shorter than the characteristic dynamical time-scale 
of the system. In that case, the effective star formation 
rate depends only weakly on the exponent n and the effi- 
ciency function /(T), and the star formation effectively 
follows the Schmidt law with the 'true' star formation 
rate depending quadratically on the gas density. Addi- 
tionally, due to the self-regulation inherent in this star 
formation prescription, the effective star formation rate 
is not strongly affected by substantial variations of / or 
C„. Therefore, we can restrict our model to the simpler 
expression 



* = C2g^ 



(5) 



We will use this expression in our analysis. 

For the description of the spatial dynamics of the 
multi-component disk we will use a ffuid dynamical ap- 
proach. In this approach, gas, stars and remnants are con- 
sidered as three ffuids, coupled by nonlinear interchange 
processes and by the common gravity. While the appli- 
cation of ffuid dynamics for the description of the gas is 



quite natural, it is not obvious that such an approach can 
be used for the coUisionless components of the disk repre- 
sented by stars and remnants. Kikuchi et al. (1997) ana- 
lyzed this question and found that the stability properties 
of disks obtained in ffuid approximation are in good quali- 
tative, and to some extent in quantitative agreement with 
the stability properties of the coUisionless models. We will 
use therefore a ffuid approximation in the analysis of per- 
turbations of a multi-component disk. The behavior of our 
model disk is described by the continuity equations, writ- 
ten for each component, a set of momentum equations, 
and the Poisson equation. In cylindrical coordinates they 
are: 



DsCTs 



Dt 



^ = (i-C)c..^ + (i-.)^ 



(6) 
(7) 
(8) 



where Dg^s,rlDt are the corresponding substantial time 
derivatives written in cylindrical coordinates: 



Dg,s.r _d ^Id ld_ 

Dt ~ dt^ rd/'''^''''^ rd6"' 



(9) 



Using equations (q) - (|8|) and deffnition (g) the momentum 
equations can be written as 

(10) 



(7s 

?7— Vs 



i>Vs 

' Dt 



'~Dt 



VPs +crsV $ + $ 



(11) 



VP^ + cr,.V $ + $ 



'H 



(l-C)C2aX + (l-^/)f V, 



(12) 



Here the surface densities, and Ug^s,r and Vg^s,r 

are the radial and azimuthal components of the velocities 
Vg,s,r of gas, stars and remnants in the disk. 

Without the exchange processes, i.e. when the right- 
hand sides of the equations ( p^ ) - (|l^) are zero, the dy- 
namics of these quantities is determined by the partial 
"pressures" of the components Pg,s,r, the self-gravity of 
the disk <& and the external gravity of the halo and bulge, 
and ^B- Mass transformations between the compo- 
nents give an additional factor for the momentum balance 
and have to be taken explicitly into account in the numer- 
ical simulations. 
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The gravitational potential <& is determined by the 
overall density of all components, and can be written in 
the form of a Poisson integral as 



K(r', 0') + as{r', cj,') + ar{r' ,cp')y dr' 



y^r^ + r'^ — 2rr' cos((/) — (f>') 



(13) 



The equation of state closes the system of equations (^ 
- (|l^). Throughout the simulations we use a polytropic 
equation of state applied to all three components: 

Pg.s.r = Kg^s,rCj],Xr (14) 

Equations are used for the analysis described in 

the subsequent sections. 

In our simulations we will use a "galactic" system of 
units in which R = 20kpc, M = IO^^A/q and the gravi- 
tational constant G is unity. With this choice we have a 
unit of time equal to 1.33 x 10^ yr, and a unit of velocity 
equal to 147kms~^. 

3. Axisymmetric Distributions 

We assume that all three components have axisymmetric 
flat rotation curves in the outer regions of the disk. The 
quasi-stationary equilibrium rotation is jointly supported 
by the gravity of the halo and bulge, the self-gravity of 
the disk and the pressure gradient: 

rnl{r) = -^{^H + ^B + <^n. + Pm) (15) 



Here Vli{r) is the angular velocity in the disk, <i>oi is the 
axisymmetric self-gravitating potential, and Pqi are the 
partial pressures of the components corresponding to their 
unperturbed density distributions. 

In our simulations, we consider two types of the equi- 
librium rotation curves and surface density distributions. 
In the first experiment, we choose the exponentially de- 
creasing surface density of the disk with density distribu- 
tions of all three components given by the expression 



(r) = Mg^s^r ■ exp - di • ^^2+^3- r'^^ ■ 



1 — exp 



( 



(16) 



with the normalization constants Mg^s,r which are the 
masses of the gaseous, stellar and remnant component. 
The halo-bulge potential determining the rotational curve 
( I^) i s similar to that used by Vauterin and Dejonghe 
( |l996D . Namely we assume that the external potential has 
the form 



H — X — = 



dr 



dr 



Mf 



(r2 



(17) 



In a thin disk with sharp boundaries the potential di- 
verges at its edges. To study to what extent the sharp 
cut-off in an unperturbed density distribution influences 
the results, we performed, additionally, experiments us- 
ing a Gaussian-type density distribution vanishing at both 
boundaries: 



{r - i?out 



1 — exp ^ — 

2s -,4 



1 — exp 



) exp 



{r - Rof 



(18) 



In this case the disk is kept in centrifugal equilibrium by 
an alternative bulge/halo distribution resulting in a radial 
acceleration of the form 



H 



dr 



dr 



B 2 



R]j 



{r^ + R%)^ 



(19) 



Both halo and bulge potentials, described by the equations 
( p7[ ) and ( pj| ) dictate the " flat" rotation in the outer region 
of the disk. The parameters Rh and Rb in equations ( [l7| ) 
and (19) determine the spatial scales of the halo and bulge 
density distributions, Voo gives the asymptotic value of the 
rotational velocity, and Mh and Mb determine the masses 
of the halo and the central bulge. 

4. Parameters of the Model 

We consider the dynamics of the disk by setting the ra- 
dius of the innermost boundary cell to i?in = 0.1, and the 
radius of the outer boundary cell to i?out = 1- 

The group of parameters r, ^, r] and C2 in the right- 
hand side of the continuity equations (||)-(||) govern the 
interchange processes between the components. Following 
Koppen et al. ( 1995 ) we choose the mean stellar lifetime 
T = lOMyr, or in our units r = 0.075. The mass fraction 
C of the newly formed massive stars was set to 0.12. This 
value corresponds to a Salpeter-IMF ranging from 0.1 Mq 
to 100 A/q and a lower mass limit of massive stars of 10 
Mq. The fraction 77 of mass ejected by massive stars back 
to the interstellar medium was taken to be 0.9. 

The parameter C2 was set to 0.1. With this choice, the 
maximum star-formation rate in our model is 0.025, or 
w ISMq/ji if the initial mass of the gaseous disk is equal 
to 0.5 in our units. This value of the star-formation rate 
is close to the maximum star-formation rate obtained in 
chemo-dynamical models for the evolution of disk galaxies 
(Samland |199|). 

We have assumed that the gas component of the disk 
is mainly composed of mono-atomic hydrogen with a vol- 
ume polytropic index 7^ = 1.67. The polytropic constant 
for the collisionless stellar component and the remnants 
was set to 2.0. There are a few arguments in favor of this 



choice. Marochnik (196(:) found that in a rigidly rotating 
disk the dynamics of perturbations can be described by 
introducing the polytropic equation of state with — 2. 
This value is consistent with the empirical "square root 
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law" found by Bottema (1993) in his studies of nearby 
spiral galaxies. He found, that the surface density distri- 
bution of stars, and their radial velocity dispersion are re- 
lated as Cs (X y^a{r). It is easy to see, that such a "square 
root law" requires the value of the effective polytropic in- 
dex to be 7s = 2. Kikuchi et al. ( 1997| ) made a detailed 
comparison of the linear stability properties of the exact 
collisionless models investigated by Vauterin & Dejonghe 
(1996) with the stability properties of this model studied 
in a fluid dynamical approach. They found a full quali- 
tative agreement between these two approaches. Thus, a 
fluid dynamical approximation can be used for the analy- 
sis of the multi-component disks. The constant Kg was set 
to be 0.04 resulting in a Toomre-stable disk (for details, 
see Sect. The values and Kr have been set to twice 
this value. This choice corresponds to a larger "sound" 
velocity of the stars, by this mimicking as well the dy- 
namical heating of disk stars as the lack of dissipation in 
the stellar component. 

All parameters including those which are not discussed 
in this section are listed in the tables in the appendix. 



5. The Code 

For solving the multi-component hydrodynamical equa- 
tions (H) ~ ( |l2|) we use a second order Van Lee r adv ec- 
tion scheme as implemented by Stone & Norman (1992) in 



a general purpose fluid dynamics code, called ZEUS-2D. 
This code was designed for modeling astrophysical sys- 
tems in two spatial dimensions, and can be used for sim- 
ulations of a wide variety of astrophysical processes. The 
ZEUS-2D code uses sufficiently accurate hydrodynamical 
algorithms which allow to add complex physical effects in 
a self-consistent fashion. This code provides therefore a 
good basis for the implementation of the nonlinear mass 
transfer processes into the multi-phase hydrodynamics. 

The Eulerian codes with the Van Leer advection 
scheme were successfully used for the investigation of the 
stability of self-gravitating disks (Laughlin & Rozyczka 
19961 , Laughlin et al. |1997|, |1998D. The main difference 



between the "standard" ZEUS-type codes and our one 
is the introduction of mass and momentum interchange 
processes between different components. These processes 
can be computed at the first sub-step of the ZEUS-type 
code which makes a generalization of the ZEUS-type code 
straightforward . 

Briefly, the code solves the hydrodynamical equations 
using equally spaced azimuthal zones and logarithmically 
spaced radial zones. For the simulations discussed here 
we mainly employed a grid with 256 x 256 cells. To ad- 
vance the solutions due to interchange processes given by 
the right-hand sides of the equations (^)-(^) we used the 
fifth order Cash-Karp Runge-Kutta routine with the time 
step limitation imposed by the Courant-Friedrichs-Levy 
criterion and the values of the parameters t and C2 in 
mass and momentum interchange processes. The Poisson 



equation is solved by applying the 2D Fourier convolution 
theorem in polar coordinates. 

6. Global Instability of a One-Component Disk 

In order to understand the influence of mass and momen- 
tum exchange processes on the evolution of global modes 
in a self-gravitating disk we performed simulations of the 
dynamics of a one-component stellar disk. The equilib- 
rium stellar disk was chosen to have the surface density 
distribution (|l6|). The parameter Ms for the stellar disk 
was selected so that the disk's total mass is equal to 0.5. 
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Fig. 1. Radial dependence of the angular velocity fi, the 
Toomre-parameter Q, and the equilibrium density a for a 
purely stellar disk (7 — 2) with exponential surface den- 
sity distribution. 



Global Amplitudes LoglO(Am) 




Fig. 2. Temporal evolution of the global amplitudes 
log(^m) (m = 1,2,3,4) for a purely stellar disk (7 = 2) 
with exponential surface density distribution. 

Figure ^ shows the equilibrium properties of the stel- 
lar disk used in the analysis. The Toomre Q-parameter, 
which is defined for the stellar disk as Q = Csk/3.36G(t, 
has a profile typical for the density distribution given by 
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Fig. 3. Contour maps of the radial velocity (in kms""'^) of a purely stellar disk (7 = 2) with an exponential surface 
density distribution at different times: t = (upper left), t = 5 (upper right), t = 10,15,20,25. The contours give 
30%, 50%, 70% and 90% of the maximum velocity in each diagram. The dotted lines correspond to negative velocities, 
whereas the solid lines give positive velocities. The zero-velocity contour is shown with a dashed line. 
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equation (p^. The Q-profile rises towards the boundaries 
of the disk, with a minimum value of 1.64 at radius 0.73 
indicating a globally stable disk with respect to Toomre's 
stability criterion Q > 1. 

The one-component stellar disk was perturbed with an 
m— armed perturbation of the form 



as{r,<j>) = cr^(r) • [1 + 0.001 cos(m<?!))] 



(20) 



1 — exp 



1 — exp I 



{r - Routf 



We chose m = 2 and m = 3 perturbations in studying the 
stability properties of our models. 

Figure || plots the development of the global ampli- 
tudes for m = 1, 2, 3, and m = 4 spiral modes in a stellar 
disk seeded by an m = 2 perturbation of the form (20). 
The global amplitudes defined by the expression 



A, 



Rou. 



a(r, (j))dre' 



-ira(f) ^ 



(21) 



illustrate the overall dynamics of the particular global 
mode. 

Figure || shows a slow exponential growth of the m = 2 
mode. The m = 1, m = 3 and to = 4 armed spirals 
grow, too, but due to the initial conditions, the to = 2 
global mode outstrips the other competitor modes during 
the whole computation. However, even at the late stages 
of evolution the amplitude of spiral perturbations is less 
then half percent, and the spiral pattern does not emerge 
from the background. The slow development of the per- 
turbation is best seen on the sequence of the snapshots 
shown in Fig. [s] illustrating the contour plots of the per- 
turbed radial velocity. (The orbital periods at the inner 
and outer boundary are 0.59 and 4.16, respectively.) 

Similar behavior was observed in the stellar disk with 
the Gaussian surface density distribution (|8|) and the ro- 
tation curve (^|) . Fig. ^ shows the equilibrium properties 
of this disk which is rather stable with a minimum value 
of Toomre's Q-parameter equal to 1.78 at a disk radius of 
R = 0.5. 

This disk was seeded with a three-armed perturbation 
of the form (pO|). As it is seen from Fig. ||, the disk develops 
a set of slowly growing modes. The behavior of the global 
amplitudes is somewhat similar to the previous case. (The 
orbital periods at the inner and outer boundary amount 
here to 0.97 and 3.65, respectively.) The m = 3 spiral 
mode has higher amplitude compared to the other com- 
petitors, and has a tendency to be saturated at the level 
\og{A-i) « —2. The perturbation, however, does not have 
any properties of a regular spiral pattern (Fig. ^ . Even at 
the late stages of the evolution the contours remain quite 
patchy, and do not resemble the regular spiral pattern. 




Fig. 4. Radial dependence of the angular velocity fi, the 
Toomre- parameter Q, and the equilibrium density a for 
a purely stellar disk (7 — 2) with the Gaussian surface 
density distribution. 



Global Amplitudes LoglO(Am) 



m=1 
m=4 




Fig. 5. Temporal evolution of the global amplitudes 
log(^m) (to = 1,2,3,4) for a purely stellar disk (7 = 2) 
with the Gaussian surface density profile. 

In the next section we will discuss how the morpho- 
logical properties of the global modes are affected by the 
change of the equation of state which is a necessary con- 
sequence of a phase transformation in a star-forming disk. 

7. Global Modes in a Multi-Component Disk 

To compare the behavior of a multi-component disk with 
the dynamics of a corresponding one-component system, 
the surface densities of all phases were initially distributed 
in accordance with equations (|l6|) or (p^). In both cases 
we choose the initial mass of the gaseous component equal 
to the mass of the one-component system. Masses of the 
"admixture", i.e. stars and remnants, were initially set to 
0.01 each which is about two percent of the initial mass 
of the gaseous component. All components were set at the 
beginning into centrifugal equilibrium, with the circular 
rotation supported by the pressure gradients of the com- 
ponents, the total gravitational field of the three compo- 
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nents and the gravity of the external halo. Fig. ^ shows the 
rotation curve and the epicycle frequency of the gaseous 
component for the initial exponential surface density dis- 
tributions. The rotation curve of the gaseous component 
is very similar to the equilibrium profiles of the stellar 
disk discussed in the previous section (Fig. ^). However, 
the gaseous component is less stable compared to the 
purely stellar disk, and the broad trough of Toomre's Q- 
parameter lies below the corresponding Q-distribution of 
the stellar disk. 




Fig. 7. Radial dependence of the angular velocity il, 
density distribution and the Toomre-parameter Q for the 
gaseous component in a multi-component disk with the 
exponential surface density distribution. 



Despite the centrifugal balance, the components are 
not in equilibrium. Mass transformations given by the 
right-hand-sides of the continuity equations (||)-(||) change 
the densities and total masses of the components, and the 
system evolves even without initial perturbations. Fig. ^ 
illustrates such mass transformation. At the beginning, 
the mass of the system is contained in the gaseous phase, 
and by the end of the computation about 90% of the gas 
has been converted into stellar remnants or long-lived low 
mass stars. The fraction of mass contained in massive stars 
drops from the initial value to 0.1% at the end of the sim- 
ulation. 

With the m — 2 perturbation given by equation (po|), 
all three phases develop a two-armed spiral pattern. Fig. 
^ shows the time dependence of the global amplitudes for 
the m — 1,2,3 and 4-armed global modes growing in the 
stellar component. Again, the m = 2 global mode prevails 
over its competitors, and compared to the purely stellar 
disk (Fig. ^ it grows about two orders of magnitude faster 
developing a nonlinear spiral pattern. 

Figures ^ and |l^ show the time sequence of the con- 
tour plots of the density distributions of the stellar and 
gaseous components of a disk with an exponential density 
profile. Initially, the density distribution evolves similarly 
to the dynamics of the one-component disk, but the sub- 









Gas + 
Stars X 



























Fig. 8. Temporal evolution of the masses of gaseous and 
the stellar component in a multi-component disk with an 
exponential density distribution. 



sequent behavior is different. A comparison of Fig. 10 with 
the simulations in a one-component disk clearly shows that 
spirals are better developed in a multi-component disk. 



Global Amplitudes LoglO(Am), Stars 




Fig. 9. Temporal evolution of the global amplitudes 
log(^,„) {m — 1, 2, 3, 4) for stellar component of the multi- 
phase disk with an exponential density distribution. 

A similar behavior was observed for the disk with the 
Gaussian surface density profile which was seeded by the 
m — 3 perturbation. Fig. shows the growth of the global 
amplitude in the stellar component of the multi-phase disk 
accompanied by the mass transformations illustrated on 
Fig. |l^. Again, a comparison with Fig. || demonstrates, 
that the growth rate of the m — 3 spiral perturbation of 
the multi-phase disk is an order of magnitude larger than 
that of the one-component system. 

The exponential growth of the m = 3 perturbation is 
changed by the lingering saturation phase occurring at the 
amplitude level log(A3) « —0.1. This nonlinear saturation 
of global modes is kno wn fo r one-componen t simulation s 
(Laughlin & Rozyczka |l996i Laughlin et al. |l997| , |l998|) . 
Our result shows that the nonlinear saturation of expo- 
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nentially growing global modes is a common phenomenon 
which occurs as well in the multi-phase gravitating disks 
experiencing phase transitions. 



Global Amplitudes LoglO(Am) 



density perturbation in remnants is more organized, and 
clearly depicts a three-armed spiral. 

The destabilizing role of the cold gaseous component 
was studied for the linear regime by various authors. Lo- 
cal analysis performed by Lin and Shu (1966), Lynden-Bell 
(| 1961| , Miller et al. ( |1970D , Q uirk (|1971| ), Jog & Solomon 
( 1984 ), Sellwood & Carlberg ( 1984D , and s emi-a nalytical 
global modal analysis of Bertin & Romeo (1988) demon- 
strated that a small admixture of gas in a stellar self- 
gravitating disk may considerably destabilize the system. 
Our simulations are in agreement with this conclusion, 
and illustrate, how the spiral structure behaves on a non- 
linear stage. Spiral structure remains well developed and 
is self-sustained after a rapid star-formation process in a 
gaseous disk when most of the gas is transformed into 
a "remnants" phase with higher velocity dispersion, and 
with a more rigid equation of state. 



Fig. 12. Temporal evolution of the global amplitudes 
log(v4m) of the m — 2 and rn = 3 modes for the stel- 
lar component of the multi-phase disk with the Gaussian 
density profile. 
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Fig. 13. Mass transformations between the gaseous and 
stellar components in the multi-phase disk with the Gaus- 
sian density distribution. 



The three-armed nature of the perturbations in the 
multi-component disk is perfectly illustrated in Figs. ^ 
and ^ which show the contour map of the surface density 
and the radial velocity of the gaseous phase in a multi- 
component disk. We note that at late phases of the evolu- 
tion the gas contains about 8% of the total disk mass, but 
nevertheless it is still a good tracer of the spiral structure. 
Moreover, our simulations allow to conclude, that the gas 
"helps" to develop spiral arms. A comparison of the den- 
sity distribution of the coUisionless remnants phase in a 
multi-component disk (Fig. |l6|) with perturbations in a 
purely stellar disk (Fig. ||) supports this conclusion: The 



8. Conclusions 

In this paper we modeled the generation of spiral structure 
in a multi-phase, star-forming disk. The gaseous and stel- 
lar phases interact as well by their common gravity as by 
mutual phase transitions due to star formation and stellar 
death. The main results of our paper can be summarized 
as follows: 

1. The multi-phase disks which undergo a gas-star 
phase transformations are unstable with respect to non- 
axial perturbations. The global spiral modes grow expo- 
nentially and saturate in a way, similar to that found in 
the one-component case. 

2. The spiral pattern grows faster and saturates on 
a higher level compared to a one-component stellar disk 
with the same mass and rotation. The spiral mode remains 
well developed, and keeps its properties unaltered if most 
of the gas is transformed into the " stellar" component with 
a stiffer equation of state. The cold gas phase remains 
a good tracer of the spiral structure even if the system 
contains about a few percent of the total mass of the disk. 

These results demonstrate that the destabilizing role 
of the cold gas component previously known for two-com- 
ponent systems with fixed background properties can be 
extended to multi-component star-forming disks which are 
allowed for rapidly changing background properties. 

Further research, however, should be performed with 
respect to self- gravitating multi-component disks. One ob- 
vious generalization is necessitated by the simplified model 
of the interaction chain used in our paper. Interactions be- 
tween the different phases in a star-forming region are reg- 
ulated mainly by the time dependent balance between the 
three components: gas, molecular clouds and stars (e.g. 
Theis et al. ( |1992D , Shore & Ferrini |1995D . Therefore, in a 
next step the gaseous phase should be split into clouds and 
an inter-cloud medium match ing th e chemo-dynamical ap- 
proach of e.g. Samland et al. ( 1997 ). Another aspect is the 
stochasticity provided by the star-formation in individual 
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clouds. Thus, in extension of Gerola & Seiden's (1978) 
analysis based on cellular automata, a study of the in- 
terplay between large-scale structure formation (including 
self-gravity) and stochastic processes would be interesting. 
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Appendix A: Parameters for Numerical 
Simulations 



Table A.l. Parameters for the exponential disk 
(cf. Eqs. (|l|) and (0)) 



General parameters: 






Courant factor 


Cf 


= 0.6 


Bulge mass 


Mb 


= 0.4 


Bulge scale length 


Rb 


= 0.4 


Halo mass 


Mh 


= 1.6 


Halo scale length 


Rh 


= 0.32 


Inner disk radius 


Rin 


= 0.1 


Outer disk radius 


Rout — 1.0 


Mass exchange parameters: 






Stellar lifetime 


r = 


0.075 


Fraction of massive stars 


c = 


0.12 


Fraction of mass returned to ISM 


V = 


0.9 


SFR parameter 


Ca -- 


= 0.1 


Disk (general): 








di = 


= 1.1 




d2 = 


= 8.0 




da = 


= 16.0 




rs = 


= 0.316 


Gaseous disk: 






Initial mass 


M, 


= 0.22 


Polytropic index 


la = 


= 1.67 


Popytropic coefficient 


Ko 


= 0.005 


Stellar disk: 






Initial mass 


Ms 


= 0.001 


Polytropic index 


"fs = 


= 2.0 


Popytropic coefficient 


Ks 


= 0.01 


Remnants disk: 






Initial mass 


Mr 


= 0.001 


Polytropic index 


Jr = 


= 2.0 


Popytropic coefficient 


Kr 


= 0.01 
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Table A. 2. Parameters for the Gaussian disk 



(cf. Eqs. (18) and (hj 



General parameters: 






Courant factor 


Cf-- 


= 0.5 


Bulge mass 


Mb 


= 0.15 


Bulge scale length 


Rb 


= 0.2 


Halo mass 


Mh 


= 10.0 


Halo scale length 


Rh 


= 0.25 


Velocity parameter 


Vao -- 


= 1.5 


Inner disk radius 


Rin 


= 0.1 


Outer disk radius 


Rout 


= 1.0 


Mass exchange parameters: the same 






Disk (general): 








= 


0.316 




Ro - 


= 0.45 




LU — 


0.05 


Gaseous disk: 






Initial mass 


M, 


= 0.5 


Polytropic index 




1.67 


Popytropic coefficient 




= 0.04 


Stellar disk: 






Initial mass 


Ms 


= 0.01 


Polytropic index 


Is = 


2.0 


Popytropic coefficient 


Ks -- 


= 0.08 


Remnants disk: 






Initial mass 


Mr 


= 0.01 


Polytropic index 


7r = 


2.0 


Popytropic coefficient 


Kr -- 


= 0.08 
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Fig. 6. Contour maps of the logarithmic surface density (Mq/pc^) for a purely stellar Gaussian disk (7 = 2) at 
different times: t = (upper left), t = 5 (upper right), 
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Fig. 10. Contour maps of the logarithmic surface density (Mq/pc^) of the stars for the multi-phase exponential disk 
at different times: t = (upper left), t = 5 (upper right), t = 10,15,20,25. The contours are logarithmic-equally 
spaced. The equilibrium model is perturbed by an m = 2-mode. 
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Fig. 11. Contour maps of the logarithmic surface density (Mq/pc^) of the gas for the muhi-phase exponential disk at 
different times: t = (upper left), t = 5 (upper right), t = 10, 15, 20, 25. The contours are logarithmic-equally spaced. 
The equilibrium model is perturbed by an m = 2-mode. 
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Fig. 14. Contour maps of the logarithmic surface density (Mq/pc^) of the gas for the multi-phase Gaussian disk at 
different times: t = (upper left), t = 5 (upper right), t = 10, 15, 20, 25. The contours are logarithmic-equally spaced. 
The equilibrium model is perturbed by an m = 3-mode. 
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Fig. 15. Contour maps of the radial velocity (in kms of the gaseous component in the Gaussian disk at different 
times: t = (upper left), t = 5 (upper right), t = 10,15,20,25. The contours give 30%, 50%, 70% and 90% of the 
maximum velocity in each diagram. The dotted lines correspond to negative velocities, whereas the solid lines give 
positive velocities. The zero-velocity contour is shown with a dashed line. The equilibrium model is perturbed by an 
m = 3-mode. 
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Fig. 16. Contour maps of the logarithmic surface density (Mq/pc^) of the stellar remnant component for the multi- 
phase Gaussian disk at different times: t = (upper left), t = 5 (upper right), t = 10, 15,20,25. The contours are 
logarithmic-equally spaced. The equilibrium model is perturbed by an m = 3-mode. 
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